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Abstract 

Understanding  the  complex  environment  of  the  rocket 
chamber  involves  good  knowledge  of  the  injection  phe- 
nomena. Understanding  the  injection  phenomena  allows 
the  rocket  designer  to  employ  time  and  cost  saving  mod- 
eling tools  to  design  a higher  performance  rocket  engine. 
The  rocket  engine  performance  is  highly  dependent  on 
the  injection  processes  within  the  chamber.  This  project 
looked  at  injection  processes  in  the  supercritical  regime 
of  the  injected  fluid,  cryogenic  nitrogen,  in  order  to  better 
understand  realistic  conditions  in  the  rocket  engines  of 
today.  The  investigation  considered  test  conditions  from 
4.0  to  6.0  MPa  at  two  different  injection  velocities  and 
temperatures.  For  the  RCM-l-A  and  RCM-l-B  cases,  the 
target  test  conditions  are  5 m/s  injection  velocity  and  120 
K injection  temperature.  RCM-l-A  is  at  4 MPa  and 
RCM-l-B  is  at  6 MPa.  Experimental  data  taken  by  Ra- 
man imaging  and  Shadowgraphy  were  compared  to  com- 
putational models  for  these  various  test  conditions.  The 
test  data  allows  comparisons  of  density,  length  scales  and 
jet  spreading  angles.  The  results  validate  the  computa- 
tional models  and  agree  with  classical  theory. 

Introduction 

There  is  much  interest  in  high  pressure  combustion  for 
the  production  of  high  power  energy  conversion  and 
thrust.  This  is  found  in  diesel  engines,  gas  turbines,  and, 
in  particular,  rocket  engines.  These  high  pressures  often 
exceed  the  critical  pressures  of  the  injected  fuel  and/or 
oxidizer.  Understanding  the  complex  environment  of  the 
rocket  chamber  in  order  to  get  the  most  power  out  of  it 
requires  a good  understanding  of  the  injection  phenom- 
ena. This  understanding  allows  the  rocket  designer  to 
employ  time  and  cost  saving  modeling  tools  to  better 
design  a higher  performing  rocket  engine.  Much  work  has 
been  accomplished  in  this  area  as  can  be  seen  in  refer- 
ences '■ 2' 3 and  4,  but  a strong  understanding  of  the  injec- 
tion process  and  the  development  of  reliable  modeling 
tools  still  requires  much  work.  This  work  takes  a step  in 
this  direction  by  looking  at  cryogenic  injection  at  pres- 
sures above  the  critical  pressure  and  modeling  of  this 
flow5.  A more  complete  discussion  of  this  effort  including 
further  testing  conditions  can  be  found  in  reference  6.  The 
comparison  between  measured  and  calculated  values 
provides  some  insight  as  to  the  reliability  of  the  modeling 
effort  and  behavior  above  critical  pressures. 

Problem 

This  investigation  looks  at  the  behavior  and  properties  of 
a cryogenic,  axisymmetric  jet  in  a supercritical  environ- 
ment. The  experiment  injects  liquid  nitrogen  above  the 


critical  point  (3.39  MPa)  through  a single  injector  into  a 
chamber  filled  with  ambient  temperature  gaseous  nitro- 
gen. The  experimental  conditions  look  at  the  effects  of 
pressure,  temperature  and  injection  speed  of  the  liquid 
nitrogen  into  the  chamber.  The  targeted  testing  conditions 
are  injection  values  of  5 m/s  and  120  K and  RCM-l-A,  4 
MPa;  RCM-l-B,  6 MPa.  Actual  injection  conditions  are 
determined  from  previous  experiments  to  determine  tem- 
perature variation  along  the  injector  tube  and  calculated 
in  the  computational  models  based  on  these  experiments. 
At  these  conditions,  Raman  and  Shadowgraph  images 
were  taken  encompassing  the  area  from  the  injector  to  60 
mm  (approximately  30  injector  diameters)  from  the  in- 
jector. The  pictures  were  used  to  determine  density  distri- 
butions, jet  spreading  angles  and  length  scales.  These 
measured  values  were  compared  to  computational  models 
for  each  of  the  cases. 

Cryogenic  Jet 

A jet  flow  has  three  distinct  zones:  potential  core,  devel- 
opment or  transition  region,  and  a similarity  region  as 
seen  in  Figure  1.  The  potential  core  contains  some  por- 
tion consisting  of  only  injected  fluid  and  reduces  in  thick- 
ness as  the  jet  mixes  with  entrained  fluid  from  the  sur- 
rounding environment. 

Self  Similar 


Figure  1 : Jet  mixing  flow  field 
After  the  potential  core  region,  the  jet  is  in  a transitional 
state  which  is  considered  the  region  of  turbulent  mixing 
for  a jet.  In  this  region,  the  energy  dissipation  and  the  jet 
behavior  tend  to  be  of  the  highest  interest  for  mixing 
purposes.  Researchers  have  indicated  the  most  significant 
influences  on  jet  development  include  the  velocity  ratio 
between  initial  jet  velocity  and  the  surrounding  environ- 
ment (Uq/iu)  s and  the  density  ratio  (p,/p„).  These  pa- 
rameters show  how  the  momentum  and  thermal  energy 
dissipates  from  the  jet  into  the  flow  field. 

At  some  distance  from  the  injection  plane,  the  jet  be- 
comes self  similar.  This  means  a function  of  only  one 
variable  can  express  the  flow  field  profiles  as  no  longer 
varying  in  the  axial  direction.  Schetz7  stated  this  to  occur 
at  approximately  x/d  > 40,  while  others  have  indicated  for 


similar  jets,  velocity  profiles  exhibit  self  similar  behavior 
as  close  as  x/d  > 20  (Schlichting*).  Other  parameters  such 
as  turbulence  intensities  (u\  v\  w’)  may  not  show  this 
behavior  until  well  after  x/d  > 200.  Also,  the  particular 
variable  employed  to  show  self  similar  behavior  varies 
between  researchers. 

Incompressible  Jet 

At  pressures  above  the  critical  pressure  and  near  critical 
temperatures,  changes  in  the  temperature  can  correspond 
to  extreme  density  gradients.  For  this  reason,  real  gas 
effects  must  be  taken  into  account  when  determining  flow 
properties.  Since  pressure  is  relatively  constant,  this  den- 
sity gradient  is  only  dependent  on  temperature.  The  ex- 
periment considers  the  density  to  be  incompressible  or 
weakly  compressible  according  to  the  classical  definition 
of  compressibility.  By  using  a real  gas  relationship  for 
density,  coefficient  of  specific  heat  at  constant  pressure 
and  viscosity,  the  computational  model  captures  the  ef- 
fects of  being  weakly  compressible  when  employing  an 
incompressible  solution  technique.  This  relationship  also 
defines  a very  strong  relationship  between  temperature 
and  density.  The  expected  results  for  temperature  and 
density  in  radial  and  axial  profiles  in  the  flow  field  are 
related  but  the  relationship  is  not  linear. 

Turbulent  Prandtl  Number  Considerations 

Schetz7  showed  the  value  of  the  Prandtl  number  actually 
varies  with  radial  location  but  on  average  is  0.7  for  simi- 
lar axisymmetric  submerged  jets.  The  evidence  from 
Pabst9  and  Sakipov"1  calculated  values  ranging  from  0.4 
to  1 .7  for  various  fluids  and  locations.  Also,  most  mod- 
eling codes  use  values  of  0.89  or  0.9  as  the  default  value 
(Wilcox", CFD  ACE  manual12).  Values  less  than  one 
correspond  to  thermal  energy  dissipating  in  the  flow 
quicker  than  momentum  energy  due  to  turbulent  mixing. 
The  resulting  non-dimensional  temperature  profiles 
would  then  be  wider  in  comparison  to  the  non- 
dimensional  velocity  profiles.  According  to  White13,  the 
value  should  be  greater  than  0.7  and  suggests  0.9  or  1.0. 
Although  this  parameter  varies  over  the  radial  profile  of 
the  flow,  it  is  common  to  choose  a constant  value  for  the 
entire  flow  field.  The  solutions  for  RCM-l-A  employed 
Prandtl  number  values  of  1.0  (White)  and  0.7  (Schetz)  as 
defined  by  the  following  equation. 
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A constant  value  provided  good  results  in  these  computa- 
tional models.  The  results  are  further  discussed  later,  but 
PrT  = 1.0  proved  to  be  more  consistent  with  our  experi- 
mental density  data  and  was  used.  It  also  allows  the  re- 
sults to  show  thermal  influences  due  to  the  changes  in  the 
coefficient  of  thermal  heat  transfer  under  these  tempera- 
ture and  pressure  conditions  without  large  influences 
caused  by  numerically  induced  turbulent  energy  transfer. 
Property  Value  Comparisons 

Non-dimensionalizing  the  flow  properties  for  the  radial 
profiles  makes  it  possible  to  compare  the  data  for  the 
different  testing  conditions  as  follows. 
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The  c subscript  refers  to  the  centerline  or  maximum  value 
for  the  profile  and  the  infinity  designates  the  environ- 
mental values.  By  this  method,  the  profiles  are  1 .0  at  the 
centerline  and  zero  outside  the  jet  itself.  Non- 
dimensionalizing  length  measurements  uses  the  jet  di- 
ameter (d).  Full  Width  Half  Maximum  (FWIIM)  values 
(r1/2),  and  axial  location  as  indicated. 

The  results  compare  the  property  profiles  at  axial  loca- 
tions of  x/d  = 0,  1.2,  5,  10,  15,  20,  and  25.  Measured 
values  at  x/d  = 0 are  not  possible  due  to  the  reflection  of 
laser  light  at  these  locations.  The  flow  properties  are  also 
compared  and  investigated  in  the  axial  direction.  Density 
and  velocity  are  of  particular  interest.  These  employ 
similar  relationships  to  the  radial  comparisons  with  one 
difference.  The  non-dimensionalization  of  these  parame- 
ters uses  the  injection  conditions  (p0,  u,,)  rather  than  the 
local  centerline  values  as  shown  in  the  following  equa- 
tions. 
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Jet  Divergence  Angle 

The  jet  divergence  angle  seems  to  be  one  of  the  most 
highly  considered  parameters  for  jet  flows.  It  lends  itself 
to  be  easily  measured  and  compared  with  other  results. 
Chehroudi  et.  al.14  provided  a comparison  of  many  differ- 
ent empirical  models  with  available  test  data  under  vari- 
ous conditions.  Of  particular  interest  to  this  experiment 
were  the  models  put  forth  by  Dimotakis"  and  Papamos- 
chou-Roshko16.  Dimotakis  investigated  the  entrainment  of 
mass  flow  into  the  growing  shear  layer  of  a free  jet.  He 
proposed  a vorticity  growth  rate  equation  seen  below 
depending  on  velocity  and  density  ratio  between  the  fluid 
flows.  For  these  testing  conditions,  the  velocity  ratio  is 
zero,  simplifying  the  following  equation  considerably. 
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Papamoschou  and  Roshko  proposed  a visual  thickness 
equation  for  incompressible,  variable-density  mixing 
layers  while  studying  the  turbulence  and  compressibility 
effects  in  plane  shear  layers.  This  relationship  uses  a 
convective  velocity  definition  to  relate  the  difference  in 
the  flows.  The  experimentally  determined  constant  (0. 1 7) 
allows  results  to  be  compared  with  axisymmetric  jet 
flows.  Again,  since  the  velocity  ratio  for  this  effort  is 
zero,  the  relationship  simplifies  considerably. 


8 =0.17 


'.-i' 


" 

( \ 

1 + 

p- 

l P J 

( \ 

( \ 

i/2  ~ 

i + 

li  o 

p. 

V « , , 

Various  methods  could  determine  the  spreading  angle 
from  the  computational  models.  Direct  evaluation  of  the 
edge  of  the  shear  layer  using  a 0.99  roll-off  point  for 
temperature,  density,  and  velocity  provides  a simple 


method  to  accomplish  this  task.  This  method  can  be  com- 
pared with  values  determined  using  a FWHM  approach. 
The  edge  of  the  shear  layer  is  difficult  to  determine  from 
Raman  images,  so  the  procedure  determines  the  location 
of  half  the  maximum  value.  The  procedure  is  to  multiply 
the  value  by  two  as  suggested  by  Chehroudi  et  al.14  to 
compare  with  visual  techniques  such  as  the  results  from 
the  Shadowgraph  images.  A similar  approach  for  the 
computational  models  also  calculate  FWHM  values  to 
use  as  a comparison  for  the  Raman  results.  The  Shadow- 
graph images  allow  direct  determination  of  the  angle. 
These  pictures  clearly  show  the  edge  of  the  shear  layer. 

Experimental  Setup 

Figure  2 shows  the  pressurized  chamber  with  the  injector 
used  in  the  experiments  presented  in  this  paper  along  with 
the  boundary  conditions  assumed  for  the  model.  The 
diameter  of  the  injector  is  2.2  mm  and  the  length  to  di- 
ameter ratio  is  greater  than  40  (Figure  3).  The  chamber 
can  be  pressurized  up  to  6 MPa  and  is  equipped  with  an 
electronic  heater  to  keep  the  wall  temperature  constant. 
Optical  access  to  the  chamber  is  provided  by  four  win- 
dows. Cold  nitrogen  is  injected  into  a warm  nitrogen 
environment  under  different  ambient  and  injection  condi- 
tions. The  temperature  of  the  injected  fluid  can  vary  from 
100  to  140  K,  the  injection  velocity  ranges  from  1 to  10 
m/s,  and  the  ambient  pressure  can  be  as  high  as  6 MPa. 

The  temperature  of  the  injected  fluid  is  generally  meas- 
ured at  position  1 (Tl,  Figure  3).  Since  the  test  setup 
includes  no  temperature  regulation  system,  the  injection 
temperature  is  varied  by  starting  the  injection  at  the  am- 
bient temperature  of  the  injector  and  the  piping.  During 
injection,  the  piping  and  the  injector  cool  down  while  the 
injected  fluid  heats  up.  When  the  temperature  of  the  in- 
jected fluid  reaches  its  targeted  value,  the  experiment 
records  the  Shadowgraph  or  Raman  images.  Since  the 
time  required  to  take  the  images  is  small  compared  to  the 
time  the  injector  needs  to  cool  down,  the  project  assumes 
quasi  steady  state  conditions. 
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Figure  3:  Injector 

The  Raman  and  Shadowgraph  images  taken  were  proc- 
essed to  produce  density,  length  scales  and  spreading 
angles  to  compare  to  the  modelled  test  conditions.  This 
procedure  is  reported  in  reference  6. 

Modeling 

The  model  for  this  multi-physical  problem  bases  itself  on 
a straight  forward  computational  approach.  The  flow  field 
calculations  employ  directly  the  Navier-Stokes  equations 
for  incompressible  flow.  Since  the  test  conditions  are  in 
the  supercritical  regime  for  nitrogen,  real  gas  nitrogen 
properties  are  necessary.  For  this  requirement,  the  model 
employs  the  Lee-Kessler17  and  Chung1 11  model  built  into 
the  CFD-ACE  software.  The  model  density  results  were 
compared  to  Younglove19  to  ensure  accuracy  of  the  input 
property  data.  Without  high  velocity  and  the  use  of  real 
gas  properties,  it  is  possible  to  employ  the  incompressible 
solution  scheme  and  still  take  into  account  the  variable 
density.  The  model  focuses  on  a steady  state  solution  to 
determine  average  property  distributions  for  this  injection 
experiment.  The  orientation  of  the  injector  also  allows  the 
assumption  of  negligible  body  forces.  A look  at  the  cal- 
culated Grashof,  Froud,  and  Reynolds  numbers  for  this 
experimental  range  shows  the  inertial  forces  to  be  the 
most  significant  with  buoyancy  and  viscous  forces  some- 
what less  reinforcing  this  assumption.  Therefore,  the 
classical  flow  equations  for  this  problem  are  as  follows: 
Continuity 


Momentum  (Navier-Stokes) 
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Figure  2:  Test  chamber 


K 


dT_ 

dx  , , 


V 


Modeling  these  equations  directly  is  very  difficult  and  is 
only  practical  for  simple  flow  situations.  For  this  reason, 
the  density  averaging  technique  derives  a time  averaged 
solution  by  using  an  average  value  and  fluctuation  to 
replace  actual  flow  parameter  values.  This  relationship 
produces  a workable  equation  set  called  Favre  Averaged 
Navier-Stokes  (FANS)  and  can  be  found  in  Cebec  and 
Smith'0  as  well  as  many  other  texts. 

The  FANS  introduces  the  Reynolds  stresses  to  account 
for  turbulence  in  the  flow.  While  several  methods  exist  to 
estimate  these  values,  the  k-e  model  seems  to  be  the  most 
appropriate.  This  method  has  shown  much  success  in 
similar  problems  and  reduces  calculation  times.  The 
model  calculations  used  this  high  Reynolds  approxima- 
tion in  the  flow  and  semi-empirical  calculation  techniques 
to  determine  flow  parameters  next  to  the  wall.  In  the 
boundary  layer  the  viscous  forces  are  much  greater  than 
shear  forces.  Launder  and  Spalding21.  The  temperature 
calculation  for  the  heat  transfer  is  treated  in  a similar 
manner.  The  program  discretized  the  governing  equations 
using  a third  order  accurate  scheme  to  capture  the  large 
density  gradients  and  damped  with  first  order  upwind 
discretization  to  maintain  stable  mathematical  computa- 
tions in  the  CFD-ACE  software  package12. 

The  software  package  takes  into  account  the  contribution 
of  heat  and  mass  transfer  from  the  turbulence  by  use  of  a 
turbulent  Prandtl  and  Schmidt  number.  The  heat  transfer 
module  solves  the  total  enthalpy  form  of  the  energy 
equation  as  shown. 

7 • [pVh  )=  V • (*Vf)+ 

dx. 

The  calculation  of  an  effective  conductivity  (K^,)  value 
takes  into  account  the  turbulence  effect.  The  following 
equation  shows  how  the  program  uses  a turbulent  Prandtl 
number  to  accomplish  this.  By  using  the  turbulent 
Schmidt  number,  the  program  handles  the  mass  diffusion 
in  a similar  way  to  calculate  an  effective  diffusion  coeffi- 
cient. 
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Grid 

The  computational  grid  used  for  this  problem  is  a struc- 
tured, 2-D,  axisymmetric  grid  with  just  over  100,000 
cells.  The  refinement  in  the  injector  region  is  critical  and 


Figure  4:  Grid 


To  show  the  solution  is  independent  of  this  grid,  the  same 
conditions  were  calculated  for  Case  3 using  several 
coarser  grids.  A grid  with  85,000  cells  compared  well  to 
the  solution  for  the  primary  grid  (100,000  cells).  The 
agreement  between  the  coarser  grid  and  fine  grid  is  very 
good.  The  axial  density  profile  exemplifies  this  agree- 
ment (Figure  5).  Therefore,  the  results  are  considered  the 
same  and  the  solution  is  independent  of  both  grids.  The 
finer  grid  was  used  for  all  the  test  conditions  to  ensure 
grid  independence  at  the  other  testing  conditions. 


Figure  5:  Density,  4 MPa,  5 m/s,  120  K 


Boundary  and  Initial  Conditions 


The  inlet  and  boundary  conditions  for  this  model  are  very 
important  and  extremely  sensitive  to  temperature.  The 
experiments  include  measurements  to  determine  the  in- 
fluence of  temperature  on  the  actual  temperature  and 
pressure  input  values  for  the  model.  Measured  mass  flow, 
temperature  and  pressure  determine  the  inlet  conditions. 
From  the  measured  values,  initial  inlet  velocity,  turbulent 
kinetic  energy  (k)  and  rate  of  dissipation  (e)  are  calcu- 
lated and  input  to  the  model. 
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The  Lee-Kessler-Chung  real  gas  properties  and  measured 
mass  flow  for  the  test  condition  determine  the  inlet  ve- 
locity from  the  density.  For  the  inlet  conditions,  inlet 
velocity  and  an  assumed  turbulence  intensity  (for  this 
model,  5%)  determine  k.  The  length  scale  (1)  in  the  tur- 
bulent energy  dissipation  rate  calculation  is  the  inlet  di- 
ameter (2.2  mm). 


The  outer  wall  of  the  chamber  is  isothermal  to  correspond 
to  the  character  of  the  wall  during  testing.  The  actual  test 
chamber  has  also  shown  some  heat  transfer  to  the  nitro- 
gen in  the  injector  tube.  To  account  for  this  in  the  model, 
the  temperature  determined  at  the  chamber  inlet  is  used  as 
the  injector  tube  inlet  temperature.  The  boundary  condi- 
tions used  are  then  an  adiabatic  wall  along  the  injector 
tube  and  the  end  of  the  chamber  at  the  plane  of  injection 
into  the  chamber.  This  allows  the  temperature  into  the 
chamber  to  match  experimentally  measured  values  for 
these  test  cases.  Also,  to  take  into  account  turbulence 
contributions  to  thermal  and  momentum  energy  transfer, 
the  turbulent  Prandtl  number  is  set  to  1.0.  This  value 
corresponds  to  Reynolds  contention  of  turbulent- 
momentum  flux  and  heat  flux  being  of  the  same  order  of 
magnitude  (White13). 

The  model  procedures  calculate  the  outlet  boundary  con- 
dition based  on  total  mass  flow.  The  chamber  is  long 


enough  to  consider  the  exit  to  be  completely  decoupled 
from  the  jet  flow.  Extrapolated  calculations  for  the  outlet 
did  not  take  into  account  the  mass  flow  being  entrained  in 
the  jet  from  the  chamber  and  therefore  gave  erroneous 
velocity.  By  establishing  exit  velocity  as  a fixed  value, 
the  calculations  could  account  for  this  problem  and 
smoothly  resolve  the  jet  flow  realistically.  The  initial 
conditions  for  the  calculations  are  simply  set  to  the  outlet 
boundary  conditions,  low  velocity  at  ambient  temperature 
and  chamber  pressure,  Figure  2. 

Convergence 

The  models  went  through  between  20,000  to  30,000 
iterations  for  the  various  test  cases  in  order  to  ensure 
convergence  of  the  solution.  Several  are  similar  to  previ- 
ously run  test  conditions  and  used  these  final  solutions  as 
the  initial  conditions  therefore  requiring  fewer  iterations 
to  converge.  A mass  flow  balance  performed  by  the  pro- 
gram provides  final  proof  of  convergence.  The  difference 
between  mass  inflow  and  mass  outflow  show  values  more 
than  five  orders  of  magnitude  smaller  than  total  mass 
flow  into  the  chamber.  During  the  iterations,  the  proce- 
dure required  the  models  to  stop  and  mass  flow  calcula- 
tions were  performed  to  ensure  proper  inlet  velocity  val- 
ues were  being  used.  The  inlet  velocity  was  adjusted  as 
appropriate  and  the  result  was  0.00%  to  0. 12  % variance 
between  model  calculated  mass  flow  and  measured  mass 
flow. 

Results 

The  comparison  of  the  various  effects  on  the  injection  of 
liquid  nitrogen  includes  density  and  velocity  profile  com- 
parisons. Experimental  data  (Raman  and  Shadowgraph 
images)  provide  us  with  information  on  density  profiles, 
divergence  angles  and  length  scales  to  compare  to  our 
model  results  as  well. 

Density  Profiles 

When  looking  at  the  progression  of  the  calculated  density 
profiles  from  x/d  = 1.2  to  x/d  = 25,  the  graphs  show  the 
development  of  the  flow  as  it  moves  toward  a self  similar 
solution.  Figure  6 from  RCM-l-A  exhibits  this  trend  seen 
at  these  testing  conditions.  The  profiles  closer  to  the 
injector  show  a flat  region  (p*  = 1.0)  near  the  center  line 
(r/r1/2=  0)  which  eventually  no  longer  exists  at  x/d  > 10. 
This  corresponds  to  the  potential  core.  Even  though  some 
of  the  testing  cases  show  little  potential  core  in  the  meas- 
ured data,  the  models  calculated  core  lengths  for  each 
testing  condition.  This  results  from  the  quality  of  the 
Raman  images.  The  Shadowgraph  images  show  a definite 
potential  core  in  these  regions  and  provide  some  validity 
to  the  model  results. 

The  transition  from  a liquid-like  jet  behavior  through  the 
transition  region  to  a fully  gas-like  jet  behavior  can  be 
easily  seen  in  the  profile  after  r/rI/2  = 1.0.  Following  the 
progression  of  the  jet  using  the  density  profiles  from  x/d 
= 1.2  to  x/d  = 10,  the  slope  progressively  increases.  For 
the  profile  at  x/d  = 15,  the  slope  has  again  dramatically 
decreased  and  the  profiles  after  this  axial  location  con- 
tinue to  decrease  slightly  and  converge.  The  profiles  for 
x/d  = 20  and  25  very  nearly  share  the  same  line.  The 
experimental  data  has  a considerable  amount  of  variabil- 


ity (no  error  bars  are  shown  for  clarity)  but  follows  the 
trend  of  the  computational  models. 


Figure  6:  Density;  4 MPa,  5 m/s,  120  K 

Figure  7 shows  another  way  to  easily  see  the  trend  for  the 
jet  to  move  through  the  transition  region  and  to  the  self 
similar  region  quickly.  This  figure  presents  density  pro- 
files (p*)  as  a function  of  r/x,  RCM-l-B. 


Figure  7:  Density,  6 MPa,  5 m/s,  120  K 
The  effect  of  using  axial  position  to  normalize  the  profile 
radial  position  shifts  the  profiles  very  near  the  injector  to 
the  right  and  dramatically  portrays  the  expected  potential 
core  (x/d  = 1.2).  As  the  profiles  are  plotted  at  intervals 
away  from  the  injector,  the  density  properties  show  simi- 
lar behavior  as  seen  in  Figure  6.  The  slope  for  the  entire 
profile  increases  until  x/d  = 10  and  then  begins  to  de- 
crease again.  This  is  opposite  to  the  behavior  in  Figure  6 
due  to  the  method  of  presenting  the  data  but  corresponds 
to  the  same  phenomena.  The  jet  has  a high  density  core  to 
an  axial  position  near  x/d  = 10,  it  then  goes  through  a 
turbulent  transition  and  then  begins  to  develop  into  a self 
similar  jet  flow.  Figure  7 also  shows  the  profiles  at  x/d  = 
20  and  25  do  not  collapse  to  the  same  line,  therefore  the 
relationship  suggests  the  flow  has  not  yet  reached  self 
similarity.  The  experimental  data  corresponds  with  this 
trend  again.  At  x/d  = 1.2,  the  very  sensitive  nature  of 
using  axial  position  to  present  radial  profiles  at  locations 
very  near  the  injector  causes  the  difference  between  the 
model  and  experimental  data. 

In  looking  at  the  agreement  of  the  density  determined 
from  the  Raman  images  and  the  calculated  model,  Figure 
8 shows  a representative  sample  of  a profile  at  x/d  = 20 
for  RCM-l-A.  The  experimental  deviation  was  calculated 
from  the  values  used  to  get  an  average  density  at  each 
location.  This  difference  is  a product  of  the  averaging 
technique  and  the  variability  of  the  data.  This  figure 
shows  a good  agreement  with  the  experimental  values. 


r/ri/2 


Figure  8:  Density;  4 MPa,  5 m/s,  120  K 

Velocity  Profiles 

The  velocity  profiles  in  the  radial  direction  lend  them- 
selves to  portraying  jet  development  from  the  injection 
point  to  a fully  developed  condition  better  than  density. 
With  the  flow  inertially  dominated,  the  expectation  would 
be  to  see  the  velocity  profiles  developing  rather  quickly 
into  self  similar  relationships.  Figure  9,  RCM-l-B  shows 
the  initial  velocity  to  be  a fully  developed,  turbulent  pipe 
flow  profile.  It  develops  into  a fairly  self  similar  profile 
very  quickly  by  x/d  = 20  and  25.  These  profiles  show  the 
edge  of  the  jet  to  be  near  r/r1/2  = 2.5. 
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Figure  9:  Velocity;  6 MPa,  5 m/s,  120  K 

The  velocity  profiles  of  the  flow  develop  much  faster 
than  temperature  or  density  distributions  for  these  mod- 
eled conditions.  Figure  10  from  RCM-l-A  shows  a typi- 
cal representation  of  the  velocity  profiles  as  well.  By  x/d 
= 10,  the  jet  appears  to  be  nearly  self  similar  although  the 
jet  is  not  fully  developed  as  seen  from  the  density  pro- 
files. 


r/x 

Figure  10:  Velocity;  4 MPa,  5 m/s,  120  K 

Centerline  Density 

The  centerline  density  calculations  prove  to  be  useful  in 
looking  at  several  aspects  of  the  jet.  Primarily,  this  axial 


property  relationship  shows  how  the  jet  dissipates  with 
distance  from  the  injector.  Observations  of  the  axial  den- 
sity profile  provide  insight  into  the  behavior  of  the  jet  as 
it  moves  through  the  various  stages  of  a jet  development. 
The  experimental  data  also  provides  a means  to  compare 
the  computational  results  with  actual  testing  conditions. 
In  Figure  1 1 Case  4,  the  computational  results  present  a 
correlation  with  the  experimental  data,  even  though  vari- 
ability in  experimental  data  makes  it  difficult  to  obtain 
exact  agreement.  If  the  temperature  is  below  the  peak  in 
specific  heat,  the  jet  is  very  dense  needing  more  energy  to 
increase  the  temperature.  Beyond  this  peak,  the  density 
ratio  of  ambient  gas  to  injected  fluid  is  higher  than  closer 
to  the  injector  and  the  jet  dissipates  rapidly. 


Figure  1 1 : Density;  4 MPa.  5 m/s,  120  K 
From  Figure  11,  the  expected  characteristics  of  the  po- 
tential core  are  very  obvious  with  the  density  ratio  ap- 
proximately constant  until  x/d  = 8.  At  this  point,  the  den- 
sity falls  off  fairly  quickly  but  does  not  reach  ambient 
values  until  much  further  downstream  (x/d  > 100).  These 
profiles  provide  insight  into  the  development  of  the  jet 
through  the  various  regions  when  viewed  on  a logarithmic 
scale  as  seen  from  Figure  12;  the  potential  core,  transition 
region  and  fully  developed  region. 


Figure  12:  Density;  4 Mpa,  5 m/s,  120  K 

Figure  12  shows  the  behavior  of  a jet  (RCM-l-A)  with  a 
very  dense  core,  the  transition  to  a turbulent  mixing  zone 
and  then  into  a fully  developed  region.  The  steep  slope  at 
x/d  = 10  to  approximately  x/d  = 30  shows  the  rapid  trans- 
fer of  both  momentum  and  thermal  energy.  At  this  point 
in  the  jet  flow,  the  density  gradient  is  lower  and  the  dissi- 
pation is  also  reduced  suggesting  a region  of  developed 
flow.  At  a point  past  x/d  = 150,  the  density  falls  off 
quickly  toward  the  chamber  value  suggesting  the  jet  has 
almost  completely  dissipated.  The  model  calculations 
show  this  location  in  the  chamber  to  be  the  far  edged  of  a 
recirculation  zone  established  to  transfer  mass  back  to- 


ward  the  injector.  This  mass  is  the  mass  entrained  in  the 
jet  flow  from  the  chamber. 

Liquid  Core  Length 

The  centerline  density  is  a power  function  of  the  axial 
distance  from  the  injector.  When  plotted  on  a log-log 
scale,  it  is  much  easier  to  determine  liquid  core  length 
and  see  the  various  regions  of  the  jet.  The  location  where 
the  density  begins  to  drop  off  rapidly  determines  the 
liquid  core  length  (xc/d).  This  location  is  easily  seen  in 
Figure  12.  For  these  test  campaigns,  several  other  condi- 
tions were  also  tested  and  calculated.  These  testing  con- 
ditions provide  information  at  various  density  ratios. 
Figure  13  shows  the  comparison  of  these  testing  condi- 
tons  for  values  of  core  length  obtained  from  the  calcu- 
lated models  against  the  Reynolds  number  relationship 
from  Harsha  and  the  relationship  offered  by  Chehroudi. 
The  two  lines  identify  the  range  Chehroudi  gave  for  the 
coefficient  in  his  relationship.  The  calculated  values 
agree  better  with  the  empirical  gas  relationship  suggested 
by  Harsha  but  tend  to  fall  at  or  below  the  lower  limit 
suggested  by  Chehroudi.  The  agreement  of  the  core 
length  calculations  to  a gas  jet  suggests  the  cryogenic  jet 
under  supercritical  pressure  behaves  very  similar  to  a gas 
jet. 


Figure  13:  Core  Length 

Unfortunately,  it  is  difficult  to  determine  highly  reliable 
values  for  the  core  length  from  the  Raman  and  Shadow- 
graph images.  Another  observation  of  the  relationship  for 
the  core  length  is  the  values  appeared  to  be  nearly  linear 
over  the  range  of  Reynolds  numbers  and  density  ratios  for 
these  experiments. 

Angle  of  Jet  Divergence,  a 

The  angle  of  jet  divergence  is  determined  from  different 
methods.  The  computational  models  used  two  primary 
methods  and  applied  each  to  velocity,  temperature,  and 
density.  The  first  one  uses  the  FWHM  values  and  deter- 
mines the  spreading  angle  from  these  radial  locations.  As 
suggested  by  Chehroudi  et  al14,  the  values  are  then  multi- 
plied by  2.  The  second  one  looks  at  the  99%  roll-off  point 
to  determine  the  jet  width.  The  computational  models 
allow  this  point  to  be  determined  quite  easily 

The  results  show  a considerable  difference  between  the 
50%  method  and  the  99%  values.  When  considering  the 
possibility  the  FWF1M  location  doesn’t  correspond  to  half 
of  the  jet  width,  another  value  for  the  2x  factor  can  be 
calculated.  For  velocity  and  temperature,  this  factor  is 
approximately  2.5  which  corresponds  nicely  with  the 
radial  profiles  when  using  r1/2  to  normalize  the  radial 


distance.  In  Figure  9,  the  profiles,  beginning  to  show  a 
self  similar  behavior,  converge  at  xkm  ~ 2.5  for  the  edge 
of  the  jet  (u*  = 0.0).  The  relationship  for  the  density  is  the 
same.  In  Figure  6,  the  profiles  converge  to  an  r/rl/2  = 3.0 
and  the  value  calculated  from  the  model  numbers  is  the 
same.  This  result  can  then  be  compared  to  the  difference 
in  the  Raman  and  Shadowgraph  methods  since  both  tech- 
niques result  from  density  relationships.  The  Raman  im- 
ages lend  themselves  to  determining  the  FWHM  values 
easier  than  finding  the  edge  of  the  jet,  while  the  Shadow- 
graphs show  the  edge  of  the  jet  fairly  clearly. 

From  the  Raman  data,  the  FWHM  difference  between 
centerline  and  ambient  density  determines  the  jet  spread- 
ing angle.  Close  to  the  injector  exit,  the  centerline  density 
is  quite  high  and  the  jet  is  quite  compact.  Added  to  this, 
the  temperature  difference  across  the  shear  layer  is  small 
compared  with  the  difference  across  the  entire  jet.  The 
dense  core  of  the  jet  therefore  determines  the  FWHM 
instead  of  the  shear  layer.  Therefore,  the  angle  is  small 
close  to  the  injector  when  determined  in  the  region  x/d  = 
0 to  10. 

In  the  images  for  x/d  = 10  to  20  and  20  to  30,  the  jet 
warms  up  and  dissipates.  The  centerline  density  is  of  the 
same  order  of  magnitude  as  the  density  in  the  shear  layer. 
Therefore,  the  FWHM  location  falls  within  the  shear 
layer.  Comparison  with  the  computational  results  shows  a 
fairly  good  agreement,  Figure  14.  For  the  coldest  cases, 
the  angles  based  on  the  FWHM  are  very  small  or  even 
negative  near  the  injector.  For  RCM-l-A,  the  minimum 
value  for  r)/2/d  corresponds  with  the  potential  core  length 
(xc/d  = 9.52).  The  conclusion  then  is  the  growth  rate  of 
the  shear  layer  and  the  radial  distance  of  the  injected  fluid 
moving  away  from  the  jet  does  not  correspond  to  the 
angle  measured  from  Raman  data  when  the  centerline 
density  is  significantly  higher  than  the  density  in  the  shear 
layer  and  ambient  gas  (the  region  near  the  injector). 
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Figure  14:  Density  FWHM;  4 MPa,  5 m/s,  120  K 

Looking  at  the  computational  model  results  and  compar- 
ing them  with  the  Shadowgraph  values  shows  a slightly 
under-predicted  trend.  In  the  Shadowgraph  images,  re- 
gions are  visible  where  the  fluid  is  in  turbulent  motion 
causing  local  density  gradients.  The  boundary  of  these 
regions  at  the  quiescent  ambient  gas  determines  the  jet 
spreading  angle.  A comparison  to  the  angle  from  numeri- 
cal calculations  using  the  99%  roll-off  point  of  the  veloc- 
ity shows  good  agreement. 

Figure  15  shows  how  the  models  for  various  density  ra- 
tios, Raman  images  and  Shadowgraph  images  compare 


with  the  relationships  put  forth  by  Dimotakis15  and  Pa- 
pamoschou  and  Roshko16.  Shadowgraph  data  was  also 
added  from  testing  conditions  using  an  injector  with  1/d  = 
1 1 .3  at  sub-  and  supercritical  conditions.  The  model  cal- 
culations seem  to  agree  fairly  well  with  the  Dimotakis 
relationship  in  this  density  ratio  region.  The  variability  of 
the  averaged  Raman  data  is  fairly  apparent  in  this  figure. 
Even  though  there  is  also  some  variability  in  the  Shadow- 
graph data,  a relationship  with  the  Dimotakis  curve  also 
seems  apparent. 


Data  from  earlier  experiments  is  also  included  in  Figure 
15.  The  experiments  at  low  density  ratios  (p«/p0)  show  a 
significant  influence  of  the  Reynolds  number.  The  Rey- 
nolds number  depends  on  the  increasing  of  the  injection 
velocity.  Correspondingly,  the  aerodynamic  forces  also 
and  are  sufficient  to  dissipate  the  jet.  At  increased  density 
ratios,  the  jet  dissipates  even  at  low  injection  velocity  and 
the  Reynolds  number  has  only  a weak  influence  on  the  jet 
spreading  angle. 
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Figure  15:  Jet  Divergence  Angle, 


Length  Scales 

The  Kolmogorov  and  Integral  length  scales  are: 


These  values  are  highly  dependent  on  the  eddy  viscosity 
method  used  for  these  calculations  to  model  turbulence. 
The  computational  models  calculated  these  values  over 
the  entire  profile  of  the  jet.  At  x/d  = 5 (Figure  16,  RCM- 
1-B),  minimum  values  for  the  integral  length  scale  (Lim) 
occur  at  the  edge  of  the  potential  core  flow  and  then  again 
at  the  edge  of  the  shear  layer.  The  later  location  of  r/d  = 
1.4  corresponds  to  the  edge  of  the  jet  based  on  99% 
change  for  temperature  and  velocity. 
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Figure  16:  LilU,  Lk0|  at  x/d  = 5;  6MPa,  5 m/s,  120  K 


The  Kolmogorov  length  scales  tend  to  be  slightly  higher 
in  the  potential  flow  region  and  somewhat  constant 
throughout  the  shear  layer.  Outside  the  shear  layer,  both 
length  scales  increase  dramatically  due  to  the  chamber 
temperature  but  have  little  physical  significance  to  the 
actual  jet  flow  itself. 

Figure  17  for  RCM-l-B  exemplifies  a typical  length  scale 
comparison  at  x/d  = 10.  Since  all  cases  show  core  lengths 
at  approximately  this  location,  the  expectation  would  be 
no  or  little  evidence  of  the  first  minimum  in  the  integral 
length  scale  as  seen  in  Figure  16.  However,  the  graph 
shows  some  evidence  of  a potential  core  region.  The 
reason  is  in  this  region  the  flow  is  in  a high  state  of  tran- 
sition to  a purely  gas  like  jet  flow.  The  second  minimum 
at  r/d  = 2.5  corresponds  to  the  values  calculated  for  the 
edge  of  the  jet  using  temperature  and  velocity. 

In  both  cases  (Figure  16  and  Figure  17),  the  geometric 
mean  length  scale  determined  from  the  Shadowgraphs 
shows  a remarkably  good  comparison  with  the  Kol- 
mogorov length  scale.  Since  these  scales  are  the  smallest 
expected  eddies  in  the  flow,  these  results  are  expected  for 
a visual  measuring  technique.  Measuring  the  integral 
length  scales  using  this  method  is  not  possible  since  these 
structures  also  include  all  of  the  easier  seen  smaller  length 
scales  as  well.  The  length  scales  closer  to  the  center  of  the 
flow  (r/d  = 0)  appear  to  increase  slightly.  For  the  colder 
flow  fields  (T  = 120  K),  the  use  of  the  Shadowgraph 


images  and  the  higher  density  gradients  in  this  region 
cause  the  calculation  technique  to  give  larger  values  for 
the  length  scales  near  the  centerline. 
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Figure  17:  Lint,  Lkoi  at  x/d  = 10;  4MPa,  5 m/s,  120  K 

Turbulent  Prandtl  Number 

The  comparison  of  the  results  of  two  different  turbulent 
Prandtl  numbers  for  RCM-l-A,  Figure  18  shows  the 
better  relationship  for  PrT  =1.0  when  considering  our 
experimental  data  for  density.  This  data  is  comprised  of 
three  individual  zones  (x/d  0..10,  10. .20,  and  20. .30  ap- 
proximately). While  the  first  region  suggests  a significant 
difference  between  the  model  and  the  measured  data,  the 
following  two  regions  corresponded  nicely  with  the 
model. 

Since  the  relationship  of  the  momentum  and  thermal 
energy  dissipation  defines  the  turbulent  Prandtl  number,  a 
value  of  one  would  cause  the  velocity  and  temperature 
profiles  to  be  very  similar  when  non-dimensionalized.  For 
values  of  turbulent  Prandtl  number  less  than  one  as  in 
Figure  18,  the  similarity  between  the  non-dimensional 
temperature  and  velocity  would  be  less  pronounced.  Our 
experimental  data  compares  better  with  an  axial  density 
profile  when  using  PrT  = 1.0.  This  suggests  the  turbulent 
momentum  transfer  is  as  important  in  explaining  the  state 
of  the  jet  as  turbulent  thermal  energy  transfer  as  sug- 
gested by  Reynolds  (White13). 


Pressure  Influence 

The  different  pressure  conditions  in  the  flow  (4  and  6 
MPa)  show  little  impact  on  the  behavior  of  the  cryogenic 
jet  since  these  pressure  conditions  are  constant  for  the 
chamber  and  are  all  above  the  critical  pressure  of  nitro- 
gen (3.39  MPa).  Figure  19  shows  the  axial  density  profile 
for  Case  3 and  Case  11,  120  K and  5 m/s.  The  Reynolds 
numbers  and  density  ratios  are  very  similar  meaning  the 
only  significant  difference  is  the  pressure.  The  difference 
in  temperature  at  the  injector  (T2)  causes  the  slight  varia- 
tion between  the  two  cases  (Case  3,  126.9  and  Case  11, 


128.7K).  At  values  so  near  the  critical  temperature  of 
126.2,  this  small  variation  influences  the  density  signifi- 
cantly. 


Figure  19:  Density;  5 m/s,  120  K 

The  influence  on  the  velocity  field  due  to  pressure  differ- 
ences also  proves  to  be  difficult  to  determine.  Since  the 
temperature  and  velocity  values  between  the  various  test 
conditions  vary  somewhat  (i.e.  RCM-l-A,  T = 126.9  K at 
uave  = 5.72  m/s  and  RCM-l-B,  T = 128.7  K at  uave  = 5.40 
m/s),  it  is  difficult  to  say  the  influence  is  only  due  to  the 
pressure  difference. 

Conclusions 

Understanding  the  complex  phenomena  of  a supercritical 
injection  flow  field  still  requires  a considerable  amount  of 
research,  but  this  work  provides  some  insights  into  sev- 
eral aspects  of  a rocket  injector.  This  work  accomplished 
this  by  examining  a single  injector  using  liquid  nitrogen 
above  the  critical  pressure.  The  various  testing  conditions 
considered  pressures  from  4 to  6 MPa  at  two  target  ve- 
locities (2  and  5 m/s)  and  a target  injection  temperatures 
of  120  K.  The  experimental  data  compares  well  with 
numerical  results  from  these  testing  conditions.  Agree- 
ment of  the  numerical  results  with  density,  length  scales 
and  jet  spreading  angles  obtained  from  Raman  and  Shad- 
owgraph images  quantitatively  validates  them. 

Under  these  testing  conditions,  the  velocity  and  density 
profiles  show  a distinct  trend  toward  a self  similar  jet 
flow  as  early  as  x/d  = 25.  The  centerline  density  profiles 
in  the  axial  direction  provide  considerable  insight  into  the 
jet  development  from  initial  core  length  through  to  the 
dissipation  of  the  jet.  The  core  lengths  determined  from 
the  computational  models  agree  fairly  well  with  the  Har- 
sha  gas  jet  empirical  approximation  and  appear  to  be 
fairly  constant  over  this  density  ratio  and  Reynolds  num- 
ber range.  Mass  fraction  profiles  show  the  importance  of 
the  coefficient  of  specific  heat,  specifically  in  RCM-l-A. 
When  injecting  fluid  below  the  critical  temperature  as 
well  as  above  the  critical  pressure,  as  the  jet  w'arms  up,  it 
will  reach  the  temperature  at  which  cp  reaches  a maxi- 
mum value.  The  heat  transfer  from  the  surrounding  gas  to 
the  colder  jet  at  this  temperature  has  a dramatically  higher 
transfer  rate  at  this  location  in  the  jet.  The  flow  conditions 
and  specifically  density  seem  to  be  most  affected  by  inlet 
temperature  variations  at  temperatures  so  close  to  the 
critical  point. 

The  various  methods  for  determining  the  spreading  angle 
proves  useful  for  comparison,  but  values  from  the  nu- 
merical models  determined  from  a 99%  roll  off  point  for 
velocity  agree  best  with  the  measured  Shadowgraph  val- 


ues  and  the  Dimotakis  model.  From  the  spreading  rate 
analysis,  a further  conclusion  reached  is  the  growth  rate 
of  the  shear  layer  and  the  radial  distance  injected  fluid 
moves  away  from  the  jet  has  no  significant  influence  on 
the  angle  measured  from  Raman  data  when  the  centerline 
density  is  significantly  higher  than  the  density  in  the  shear 
layer  and  the  ambient  gas.  This  result  is  apparent  from  the 
FWHM  method  of  determining  the  spreading  angle  from 
the  Raman  images.  The  resulting  angles  do  not  corre- 
spond with  the  Shadowgraph  or  the  numerical  results 
very  well.  The  procedure  of  multiplying  the  FWF1M  val- 
ues by  a factor  of  2 also  does  not  prove  to  be  very  accu- 
rate when  applied  to  the  numerical  results. 
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Nomenclature 

C,  Constants 

cp  Coefficient  of  specific  heat  constant  pressure 
d Injector  diameter  (2.2  mm) 

g Universal  gas  constant 

h,  H Enthalpy 

J Momentum 

k Turbulent  kinetic  energy 

kT  Eddy  conductivity 

k;,T  Effective  kinetic  turbulent  energy 

K Thermal  conductivity 

Lint  Integral  length  scale 

Lkol  Kolmogorov  length  scale 

PrT  Turbulent  Prandlt  number 

r Radial  location 

r„2  Radial  location  of  half  maximum  value 

of  a parameter 
T Temperature 

u Axial  velocity 

u’,  v’,  w’  Turbulence  intensity  values 
x Axial  location 

Y Mass  fraction 

a Jet  divergence  angle 

5,,,’  Vorticity  growth  rate  (radians) 

e Turbulent  energy  rate  of  dissipation,  optical 
efficiency 

k Thermal  conductivity 

v Kinematic  viscosity 

vT  Turbulent  eddy  viscosity 

p Density,  number  density 

t Shear  stress 

Sub-  and  Super-scripts 
ax  Axial  direction 

c Centerline  value  at  specified  axial  location 

i,  j,  k x,  y and  z,  Cartesian  directions 

o Centerline  value  at  the  injector 

rad  Radial  direction 

T Turbulent  parameter 


°°  Chamber  property  away  from  the  jet 
*,  + Dimensionless  values 

‘ Fluctuation  from  the  mean  value 
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